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1. Introduction 

Physical observables are gauge invariant, so it does not seem unreasonable to assume that 
it should be possible to bypass the unpleasantness of gauge fixing in any fundamental 
approach to gauge theories. Indeed, this is one of the virtues of the lattice approach and 
for important classes of "gold-plated" observables, such as masses or decay constants of 
stable hadrons, there have been spectacular successes jlj]. However, not all observables fall 
into this class and a basic fact of life is that for many important observables, being gauge 
invariant does not mean that the issues concerning gauge fixing can be avoided. Indeed, 
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important properties of some observables only become clear once their connection to gauge 
fixing has been clarified. A fundamental example of this, and one that we shall focus on 
in this paper, is colour ||. 

Although it sounds paradoxical, gauge fixing offers a route to the construction of certain 
gauge invariant quantities. In particular, charges constructed in this way || have better 
infrared properties than conventional fields and offer a possible resolution to the infrared 
crisis that plague gauge theories The main virtue of this approach to charges is that 
they are not simply made gauge invariant and then assumed to overlap with some unknown 
physical state (such as in the use of Wilson lines to probe the inter-quark potential on the 
lattice) but that their physical significance is built into their very construction from the 
beginning. 

The process of using gauge fixing to construct a gauge invariant charge is called dress- 
ing. Physically, we envisage the matter field being dressed by gauge fields to make it gauge 
invariant. This, though, is done in a very precise way which reflects the kinematics of the 
charges involved. So the dressing for, say, a static charge is different to that of a moving 
charge: each one will have a special gauge fixing that is adapted to it. For a static charge 
it is the Coulomb gauge that plays the central role and in that gauge the gauge invariant 
static charge takes its simplest form. 

Hence for observables such as colour the thorny issues associated with gauge fixing take 
on a new urgency. In particular, the Gribov copies || associated with a wide class of gauge 
fixings cannot be treated as just a technical issue to do with over-counting configurations, 
but will have direct physical consequences. In [0], general arguments were given for the 
global breakdown of colour charges in QCD due to such copies. This implies that Gribov 
copies and the underlying global topology of the Yang-Mills configuration space ||, lead 
inexorably to the absence of colour in unbroken gauge theories. Many advances are still 
needed to fill out this route to confinement. The infrared dynamics of these charges is subtle 
and although there has been some progress ||, ffl, much still needs to be done, particularly 
in understanding collinear divergences. However, recent lattice simulations Q have used 
such dressings to calculate the non confining inter-quark potential, which gives additional 
support to this programme. 

In this paper we want to expand on the arguments given in Q and understand precisely 
how Gribov copies obstruct the construction of a charge in a non-abelian gauge theory. To 
this end it is desirable to have a set of simple examples of such copies. Although the 
literature on this topic dates back to Gribov's original 1977 paper (for a recent review 
see ||), we should note that charges are never considered and that the bulk of such work, 
see e.g. |15|| , deals with gauge transformations U(x) that are non-trivial asymptotically, 
i.e., U(x) -t* 1 as \x\ — ► oo and as such correspond to 'large' gauge transformations which 
belong to the disconnected part of the group of gauge transformations. 

Due to this restriction, simply adding charges to what has been done before is not as 
simple as it sounds since we have shown Q that we must asymptotically have U(x) — > 1 
(or at least the centre) for a colour charge to be well defined. Two papers are worth noting 
in this context. Henyey Q constructed an explicit set of copies of an axially symmetric 
gauge field, while Grotowski and Schirmer [ 10 1 proved that a large class of spherically 
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symmetric configurations have copies. In both cases the gauge transformations are regular 
ones for which a colour charge can be defined. Although the Grotowski and Schirmer 
paper gave no examples it is very suggestive as it offers the possibility of a large class of 
explicit configurations that would allow us to probe the extent to which colour charges 
can be defined. It is important to note that in most discussions in the literature the 
coupling g is hidden. We will include the coupling as this will allow us to demonstrate the 
non-perturbative nature of the copies even when they are generated infinitesimally. 

The plan of this paper is as follows. After this introduction, in Section 2 we will recap 
how colour is defined and the significance of this for the allowable gauge transformations. 
We will then present the dressing approach to the construction of static colour charges 
which extends the methods given in In Section 3 we will enlarge on the arguments 
given by Grotowski and Schirmer [10], and then construct explicit examples of spherically 
symmetric copies. The role of the centre in these constructions will be also be investigated 
and this will throw some light on the complicated topology of the Gribov horizons. 

Following this, in Section 4 we will investigate the extent to which such configurations 
can have copies within the Gribov horizon. Our results here contradict one of the con- 
clusions in [1C]. We will conclude in Section 5 by tracing the impact of the copies on the 
colour charges constructed in Section 2. We will see explicitly that the colour acquires a 
gauge dependence which emerges through the copies. Although this breakdown of colour 
will be manifestly non-perturbative, we shall see that it can arise even infinitesimally. After 
a general conclusion we include several appendices where technical results are presented. 



2. A world of colour 



Imagine a world of colour. By this we mean that there exist (asymptotic) states with well 
defined colour. The fundamental question at the heart of confinement is then: to what 
extent is this possible? Charges in gauge theories are notoriously hard to construct, even 
in the abelian theory [|J|. However, we shall see that there is an unavoidable obstruction in 
unbroken non-abelian theories which stops the partonic, effective colour degrees of freedom 
from becoming genuine asymptotic states. 

To proceed let us fix notation and recall some basic facts about gauge theories. The 
gauge transformations on the vector potential and matter fields are: 

A u -> A u a := U~ l AJJ + -U- l dnU and if; -» i(j u := U^ip , (2.1) 

g 

where = A^T a and T a are anti-hermitian matrices in the appropriate representation 
of SU(iV) and g is the coupling constant. Writing U = ex.p(ga a (x)T a ), we identify the 
generator of gauge transformations with the charge 

G(a) = J d 3 x G a (x)ga a (x) := J d 3 x (j a (x) - i(A^) a (*)) <?a a (z) , (2.2) 

where D{ = d{ + <7p4j,-] is the covariant derivative, Jq = —i^oT "^ and Ef(x) is the 
chromo-electric field. The Gauss law constraint used to identify physical states is then 
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taken to be 1 

G a {x)\ physical) = 0. (2.3) 

2.1 Defining the colour charge 

The conserved Noether current associated with the gauge transformation is given by 

j°(x) = J£(s) - f abc F b v (x)A» c (x) , (2.4) 
where f abc are the structure constants. Hence the colour charge Q = Q a T a is given by 

Q a = j d 3 x f Q (x) = J d 3 x (j » - f<*°A$(x)E?(xj) . (2.5) 

This is conserved but it is not gauge invariant and hence cannot immediately be identified 
with a physical observable. However, acting on physical states defined by fl2,3| ) it takes a 
particularly simple form 

Q\ physical) = - / d 3 x diEAx)\ physical) = - lim / ds ■ E(x)\ physical) , (2.6) 

9 ./ gR-*°°Js 2 R 

where we have used the divergence theorem to get an integral over the spatial two-sphere of 
radius R. In this form it is straightforward to apply a gauge transformation to the charge 
to get 

Q u = -lim [ ds-U~ l (x)E(x)U(x). (2.7) 

9 R->oo J S 2 r 

If we require a physical colour charge then we must have Q u = Q when acting on physical 
states. To ensure this we must have that U (x) tends to a constant Uoq in the centre of the 
group. So for SU(2) we have a physical colour charge as long as 

U(x)^±t as \x\^oo. (2.8) 

2.2 Constructing colour charges 

Having clarified what the colour charge is and the states that it must act upon, we now 
need to construct physical states that carry colour. An immediate observation is that the 
matter fields ip(x) are not suitable since they are not gauge invariant and hence neither can 
be identified with observables nor have a colour. What is required is a composite operator 
of the form 

V(x) = h~ 1 (x)il}(x), (2.9) 
where the field dependent dressing h~ 1 (x) transforms as 

hr x {x) -^h~ x U{x), (2.10) 

under the gauge transformations so as to compensate for the transformation of the matter 
field. 



1 The relationship between this condition and the gauge invariance of physical states is a subtle one that 
we do not need to worry about in this paper. 
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Dressings may be generated using gauge fix- 
ing conditions and why this should be so is quite 
easy to see. Recall that gauge fixing is supposed 
to pick out a unique representative from the or- 
bit of gauge related potentials. So let x(A) = 
be a gauge fixing condition. For any configura- 
tion A there is a unique gauge transformation 
h = h[A] which takes the representative A to the 
point on the orbit where the gauge fixing con- 
dition holds, i.e. x{A h ) = 0. Now consider a 
different point on the orbit, A u , then there will 
be a dressing h u = ^[^4^] that takes it to the 
gauge fixing condition. From this it follows that 
x((A u ) hU ) = x(A UhL ) = and therefore by uniqueness that h u 
transforms as a dressing. 




Figure 1: The 'hoped for' view of gauge 
fixing - that the gauge orbits intersect 
the chosen gauge slice (e.g. diAf = in 
Coulomb gauge) once and only once. 



U 1 h so that h 



-i 



As discussed in [11], there are various possible choices for the dressing but most of 
these would represent highly excited states with little or no overlap with a physical charge. 
The dressing may be fixed by imposing the condition that the resulting charge is static. 
In that case, the dressing (or at least the part responsible for making the matter field 



gauge invariant, which is all which will concern us in this paper, see however, [13]) can be 
identified with the field dependent gauge transformation which takes a potential A^ into 
Coulomb gauge di{A h )i = 0. Tests of this construction have been carried out by studying 
the infrared properties Q of the fields and the inter-quark potential [12, between two 
such charges. 

So, for a static charge we need to solve the Coulomb gauge condition for h, 



di(A h )i = d l {h~ 1 A l h + hi- x dih) = . 

To help in this we note that the above may be written as 

d i {A l + -h- 1 D % h)=0, 
9 



(2-11) 



(2.12) 



where now the covariant derivative acts on the group elements by simple commutation: 
D i h = d i + g[A i ,h]. 

If we write h = e v then, as discussed in Appendix A, we can write 



1 — e v 
hT x D,h = DiV. 



(2.13) 



where we define powers of v acting to the right on some operator B by v°B = B, ~~vB 
[v, B] and ~v n B = ~v(~v n ~ 1 B). Acting on this by di we get 



di(h- l Dih) 



1 - e v 



(diDiv) + 



— =; — ' —°i v 

V V 



D iV . 



(2.14) 
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Substituting this expression into (|2.11 ) allows us to construct the static dressing, at least 
perturbatively. In particular, if we expand v as a power series in the coupling v = ^2g n v n , 
then we readily find that, for example, 



vi 



djAj 
V 2 



d 

and v 2 = ^([v u Aj] + ^[djV U vi]) . 



(2.15) 



In this way we can perturbatively identify the static quark with the dressed field 

■qj = h- 1 ^ = e 9vi+g 2 v 2 +... 1 p (2.I6) 

In this and what follows the Laplacian is V 2 = <9i<9j. 



3. Gribov copies 

The simple view depicted in Figure 1 of gauge fixing as selecting a unique representative 
from each orbit does not hold in practice and, as shown in Figure 2, the actual interplay 
between the gauge orbits and the gauge slice can be quite complicated. The existence 
of such copies, even when the gauge transformations are restricted to lie in the centre 
asymptotically, are guaranteed by Singer's theorem ||. Here we shall show how to explicitly 
construct such Gribov copies. 




o 



d,A] = -1 



Figure 2: The Gribov problem: a gauge field A may be gauge equivalent to multiple copies in the 
same gauge slice as the gauge orbit A u intersects (points 1,2,3. . . ) or becomes tangential to the 
gauge slice (point 4). 



3.1 Spherical copies 

In the remainder of this paper we will take our gauge group to be SU(2), where 

[T ,T \ =e T , T = —a , 
with a a the Pauli matrices. We start with configurations of the form 

A^x) = /(r) ~ 1 g iafe - . (3.1) 
r r 

This is the simplest type of spherically symmetric configuration and is the form consid- 
ered by Grotowski and Schirmer 113]. The key point to note is that this configuration 
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automatically satisfies the Coulomb gauge condition diA^ = 0. Consider now the gauge 
transformation described by 



(3.2) 



U(x) = exp \2ga{r)—T a 

x^ 

= cos (ga(r))t + isin (ga{r)) — a a , 

r 

where we are setting a a (x) = a(r)x a /r. This transformation takes A to A u where 

{A u ) a % = (f(r) cos(2 5 a(r)) - 1 + 2(1 - g' 1 ) sm 2 (ga^j (3.3) 

, x l x a sin(2gq) , _^ 
-2a (r)— (f{r)-l + g 



This is an example of a more general spherical configuration. We now impose the condition 
di(A u )® = Oto find copies. Only the final two terms in ( |3.3| ) make a nontrivial contribution 
to the Coulomb gauge condition, so we need to solve: 

+ «(/M - 1 + ,-') (f - ^)) = 0. (3.4) 

From this, using the properties of the transverse projector, we rapidly arrive at the key 
equation for the existence of copies: 

rV'(r) + 2ra(r) - (/(r) - 1 + g~ l ) sm(2ga(r)) = . (3.5) 

This is the equation that we must solve to find the initial configuration, /, and the gauge 
transformation, a, which takes us to a copy. Note that Grotowski and Schirmer proved the 
existence of solutions to this equation (at g = 1) along with suitable boundary conditions. 
As such we need to discuss the boundary conditions that we use (which are more general 
than those considered in [[To|] ), 

For (3^2) to be a regular gauge transformation we require that ga(r) — > nir as r — ► 
for neZ, We may take n to be even or odd so that U{x) approaches ±1 at the origin. 
As we have seen, to define colour we need similar boundary conditions at spatial infinity. 
Our allowable boundary conditions are therefore 

qa(r) — > < , ^ as r —> 0, oo =>■ U(x) — » < . (3.6) 

w [ (2n + l)vr K ' | -1 y J 

Boundary conditions on the function / come from finite energy and norm restrictions. The 
L 2 -norm of a configuration A is defined in the obvious way as 



Applied to our class of configurations ([O]) we see that 

POO 

\\A\\ 2 = 8vr / dr (/(r) - l) 2 . (3.8) 
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So if we require a finite norm, given that / is regular, we must have that 



f(r) — > 1 as r — > oo . (3-9) 

There is no compelling physical reason for imposing such a finite norm condition on our 
configurations. However, following Grotowski and Schirmer we shall do so here as mathe- 
matically it leads to simplifications. Note, though, that weakening this condition is possible 
and we will discuss this in Section 4.4. 

In contrast to the norm of the configuration, the norm of the field strength is of 
central importance as it is the energy of the configuration. For our spherically symmetric 
configuration we see that 



- (/ 2 (r) _ 1)2 
2 



||F A || 2 = 4tt / dr U ^ 2 ; + 2/'(r) 2 . (3.10) 



Finiteness here implies 

f(r) — * 1 as r — » , 

V ( 3 - n ) 

/ (r) — ► as r — > oo . 

Note that, as expected, the copy configuration will have the same energy but a different 
norm given by 

\\A U \\ 2 = 4vr/ dr 2 [/(r) cos(2ga(r)) - 1 + 2 sin 2 ( 5 a(r))(l - g' 1 )] 

Jo (3.12) 

+ 2sin 2 (2 5 a(r))(/(r) - 1 + cT 1 ) 2 + 4r 2 (a'(r)) 2 . 

The rate at which the boundary conditions ( |3.6| ), ( |3.9| ) and ( |3.11 ) are attained is dictated 
by the convergence of the norms of A, Fa and A u . 

3.2 Explicit examples 

To solve (3.5) we initially take the boundary conditions used in |1C] where a — > for both 
r — > and r — ► 00. The idea here is to use Henyey's trick Q of inputting the gauge 
transformation, a, by hand and then finding the copied configuration, /, via 

/(r) = r VfeI±WW +1 _l. (3 . 13) 
sin(2pa(r)J 5 

Now comes a simple observation: using the limit xj sin(x) — > 1 as x — > and given the 
boundary conditions on a, we can automatically satisfy the boundary conditions on / if 
we insist that 

r 2 a"(r) + 2ra'(r) — > 2a(r) as r — > and r — > 00 . (3-14) 

The equation 

r 2 a"(r) + 2ra'(r) = 2a(r), (3.15) 
is easy to solve resulting in two types of solution: 

a(r) oc r and a(r) oc — . (3.16) 
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Note that these relations need to be interpreted with a little care as a is dimensionless. 
More properly we should write, for example, that a(r) oc r/ro where ro is some fixed length 
scale. For simplicity, we shall omit reference to this scale. 
Hence we see that if 



a(r) 



as r 



^ as r 



0, 

oo , 



(3.17) 



then we satisfy both the boundary conditions on a(r) and those on /(r). These are the 
only conditions we must satisfy in order to find finite norm, finite energy configurations 
with Gribov copies, and there are very many ways to solve ( |3.17j ). The simplest examples 
are where a(r) is a rational function. For example, take 

a ( r ) = -J—, (3.18) 



1 + H 



The resulting f(r) derived from ( |3.13 ) clearly has a coupling dependence. We find 



f(r) 



2r 



1 - 1r A + r fc 



(1 + r 



3^3 



sin 



2gr \ 



+ 1 



1 

9 



(3.19) 



The non-perturbative nature of these solutions is now explicit in the unavoidable presence 
of factors of g in the denominators. In Figure || we plot two such configurations. 

The norms of this configuration and 
of its copy, for example at g = 1, are 

HAH ~ Vl37 ~ 11.7, 



'136 ~ 11.6. 



and the energy is \\Fa\\ — v340 ~ 18.4. 
Note that the set of connections satis- 
fying the Coulomb condition is a linear 
subspace of the Yang-Mills configuration 
space A, which is itself a convex space. 
Thus we can meaningfully use the norm 
to measure the difference between con- 
figurations. In particular, the distance 
between the copies is 




Figure 3: f(r) as constructed in ( 3.18| ) and ( 3.19| ). 
The upper curve has g — 1, the lower g = 1/2. 



47T 



dr 2 [/(r)(cos(2flia(r)) - 1) + 2sin 2 (#a(r))(l - cT 1 )] 2 
+ 2sin 2 (2 5 a(r))(/(r) - 1 + g~ l f + 4r 2 a'(r) 2 . 



(3.20) 



Now, following 14], if we introduce a parameter f3 via the replacement a{r) 
it is clear that 



||A^-A|| 



0. 



f3a(r), then 
(3.21) 



> as 13 

That is, the copies coalesce as (5 gets smaller. For j3 very small the two non-perturbative 
copies are infinitesimally close and we are clearly approaching a configuration on some sort 
of boundary, a point to which we will return in the next section. Consider the example 
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a r 



r/3 



(3.22) 



1 + r 3 ' 

Provided that |2^a(r)| < tt (we will re- 
turn to this point later, for the moment 
this is easily satisfied, especially as we 
will shortly take the limit of f3 — ► 0) there 
are no divergences in f(r) and as f3 — ► 
the gauge configuration and its copy co- 
alesce, see Figure |j. The corresponding 
f(r) in this limit is 



. 5 
. 4 

0.3 
0.2 
. 




f(r) 



1 



9r d 



with norm 11^41 



g(l +r3)2 



8tt 



(3.23) 



Figure 4: As (3 — > 0. 

a(r) is plotted for /3 = 



a(r) approaches 0. Here our 
= 1,1/2, 1/5, 1/20. 



3 3 /4g . The energy of this configuration is 



\Fa\ 



32vr 2 (10-21 5 + 2V) 



3.3 The centre 



We have so far considered only those transformations with a(0) = a(oo) = implying 
U(x) approaches unity at zero and spatial infinity. Recall, though, that in order to have 
a regular transformation we must have ga(r) —* nir as r — > 0, so that ?7(x) — ► ±1 there, 
and that we may relax the boundary condition at infinity to U(x) — > ±1 without losing 
the ability to define a physical charge. 

The gauge transformations we consider there- 
fore fall into two classes, best summarised by their 
boundary conditions, see Table 1. The examples 
considered previously are of 'closed' type - this 
language will shortly be clarified. We now wish to 
give solutions to ( |3.5| ) for general open and closed 
gauge transformations. 

We begin by expanding the discussion of the limits previously used to solve (3.5). For 
either open or closed transformations we may write ga = nix + gu{r) such that u)(0) = 0. 
The periodicity of sine in f(r) leads to the same problem as previously considered, 

r V(r) + 2ruj'(r) 





r -> 


r — > oo 


'Closed' type 


U^±t 


U -> ±1 


'Open' type 


U^±t 





Table 1: Taxonomy of open and closed 
boundary types 



sm(2gu(r)) 



+ 1 



1 

9 



where a{r) is now replaced by u(r). We conclude that uj(r) ~ r for r small. Now write 
got = rrnr + g$l{r) such that O(r) vanishes at infinity, then we must solve 



Sir) 



r 2 n"(r) + 2rn'(r) 
sin(2#0(r)) 



1 

+ 1-- 

9 
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and we see that f(r) satisfies the boundary conditions ( |3,9| ) if J7(r) obeys the same asymp- 
totic behaviour as the original a(r), that is J7(r) ~ 1/r 2 at infinity. Therefore all our gauge 
transformations must behave like 

a(r) = mr + const, r + . . . as r — > 0, a(r) = mir + const. r~ 2 + . . . as r — > oo. 

Now that we have identified this shared small r and asymptotic behaviour we will discuss 
the two classes of transformations in more depth. 

3.4 Open configurations a Q 

'Open' configurations interpolate between U(0) = ±1 and {/(oo) = =Fl. Therefore a (r) 
obeys 

f 2nvr, f(2m+ l)vr, 

goL (r) —* < as r — > and ga [r) —* < as r — > oo, (3.24) 

[(2n + l)vr [2m7r 

and clearly must pass through an integer multiple of 7r/2 where /(r) diverge unless we 
construct it carefully. The limits we employ at r = and r — > oo may also be applied here. 
Suppose ga (c) = n/2. Expand ga a (r) = ir/2 + gj(r) with 7(c) = 0, then 

rV(r) + 2rV(r) 1 
/(r) " sm(2 57 (r)) +1 5 " (3 - 25) 
Aside from an initial minus sign, our problem is unchanged from previous arguments. We 
know that this function will be well behaved at r = c if 7(7*) is a linear combination of r 
and 1/r 2 near r = c. The condition 7(c) = implies 

j{r)~x(r-^\ (3.26) 

near r = c for A some constant. A similar trick works at all multiples oin/2. This identifies 
the behaviour of a Q (r) at all potentially problematic interior points. 

Given an f(r) and the a Q (r) which generate copies we may consider the gauge trans- 
formation a Q (r/P) in the limit j3 — > 0. The distance between copies is 

poo 

\\ A U(*o) _ A f = 4n dr 2 [/( r )( cos (2 5 a (rr 1 )) - 1) + 2sin 2 ( 3 a )(r/r 1 )(l - g' 1 )) 2 
Jo 

+ 2(/(r) - 1 + g- 1 ) 2 sm 2 (2ga (r(3- 1 )) + 4r 2 (3 2 <a' (r /T X ) 2 • 
Changing variables r = v/3 yields 

/•oo 

pC/K) _ A ^2 = dv 2 [ f (p v )( C0S ( 2 ga (v)) - 1) + 2 sin 2 (ga (v)) (1 - (T 1 )] 2 

/o 

+ 2(/(/3t>) - 1 + g- 1 ) 2 sm 2 (2ga (v)) + 4i> 2 (cU (t;)) 2 . 

(3.27) 

As (5 — > 0, f{(5v) — ► 1 for all u and the integral above becomes 

/ dw 2 [cos(2 5 a (w)) - l + 2sin 2 ( 5 a (t;))(l -g -1 )] +2g~ 2 sm 2 (2ga (v))+4v 2 (d v a (v)) 2 , 
Jo 

which is seen to be finite using the behaviour ga (r) ~ it + const. jr 2 at infinity. As the 
integral is multiplied by f3 the distance between the copies vanishes as — > 0. 
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In this limit all features of a are lo- 
calised close to the origin and U(x) — > 
U(oo) in the open interval r G (0, oo), 
whereas f7(0) = — U(po) independent of 
(3. This is illustrated in Figure |fora o (0) = 
0, ga (oo) = 7r, so that i7(0) = 1 in- 
dependent of (3, but U(x) — > — 1 for all 
x / as /? -> 0. 

Compact expressions for open config- 
urations seem hard to find, however in 

Appendix B we show how they can be Figure 5: As (3 — > 0, ga {r/(3) flows to the left, 
constructed by patching functions together, approaching tt for all non zero values of r. 




3.5 Closed configurations a c 

Closed type configurations a c interpolate between U(0) = ±1 and U(oo) = ±1. When the 
a c approach the same multiple of tt at both zero and infinity they are, in general, given 
by translations of the a{r) we originally considered (for the special case a(0) = q(oo)), of 
the form ga c (r) = (2n + l)ir ± ga(r). Provided the original a(r) was suitably bounded so 
that it did not cause a divergence in f(r) then neither will a c (r), in fact these result in the 
same f(r) and copy. 

For example, including the parameter /?, an example of a closed transformation is 



ot c (r) = pa{x) = — — -. 

g g 1 + r d 



(3.28) 



This ga c approaches tt at both r = and r = oo and f(r) is well behaved for all (3 < 1. 
As (3 — * the gauge transformation approaches U(0) = U(oo) on the (half) closed interval 
[0, oo). The distance between copies is, using the translation properties of the trigonometric 
functions, unchanged, 

\\A U ^-A\\ 2 = \\A U ^ - A\\ 2 
r°° 

= 4tt dr 2 [f(r)(cos(2gpa(r)) - 1) + 2 sin 2 ( 5 /3a(r))(l - g' 1 )} (3.29) 
Jo 

+ 2sm 2 (2g(3a(r))(f(r) - 1 + g' 1 ) 2 + 4r 2 /3 2 a'(r) 2 , 

and again the distance between copies vanishes as [3 — > 0. 

If a c interpolates between different odd multiples of tt then, provided a c behaves as in 
(3.26) at multiples of vr/2, we may again consider the transformation a c (r/(3) as (3 — > 0. 
Here, as in the open case, all features of a c (r) become localised around r = and U{x) 
approaches U(oo) at all non zero points. The difference between this and the open case is 
that U(0) = U(oo) from the outset, so that although a c only approaches a single value on 
the open interval r £ (0, oo), the gauge transformation U(x) — > U(0) = U(oo) again on the 
(half) closed interval r G [0, oo). 
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4. The Gribov horizon 

To clarify what is going on here we need 
to recall some basic definitions associated 
with Gribov copies. Given a configura- 
tion with its copies, the associated Faddeev- 
Popov operator is commonly used to pick 
out a subset of copies. 

In the Coulomb gauge, with diA^ = 
and U = 1 + a a T a an infinitesimal gauge 
transformation, the equation for the exis- 
tence of copies becomes 



2.5 








d i (u- 1 A i u + -u- 1 d i u) 

9 

d i {A i -[a a T a ,A i ] + -d i a a T a + 
9 



0.5 



Figure 6: The small (3 behaviour of the closed type 
configuration ( 3.28 ) with .9 = 1. The corresponding 
gauge transformation U(x) — ► — 1 for all x € K 3 as 
/3->Q. 

1 



--V A a a T a 
9 



where T>a = —diDi is the Faddeev-Popov operator. We see that the condition for the 
existence of a copy A u close to the original configuration A is the existence of zero modes 
in the Faddeev-Popov operator. The first Gribov region is defined as the subset of config- 
uration space where the Faddeev-Popov operator is positive definite. It is bound by the 
'Gribov horizon', on which the Faddeev-Popov operator picks up its first zero mode. Out- 
side this region, see || for a fuller discussion, the zero mode becomes a negative eigenvalue 
mode and further zero modes of the operator appear. We may further decompose con- 
figuration space into Gribov regions defined by the number of negative eigenvalue modes 
of the Faddeev-Popov operator 2 . Restricting gauge fields to lie within the Gribov hori- 
zon is common practice in, for example, lattice calculations (where Coulomb gauge fixing 
is implemented via a variational argument) but does not eliminate copies, nor is there a 
physical motivation for this restriction. In this section we will consider the proof in [10] 
that spherical copies lie inside the Gribov horizon. 

4.1 The Grotowski-Schirmer argument 

We begin with an overview of the Grotowski and Schirmer argument. The Faddeev-Popov 
operator in Coulomb gauge is 

V A 



-d i D i = -V 2 -g[A i ,d i -] 



(4.1) 



although for now we restrict to g = 1, as in [10]. It should be noted that the authors of 
|1C], to which we will frequently refer, use A for both the full Laplacian and only the radial 
part in spherical co-ordinates. Given an arbitrary state we may decompose it in spherical 
vector harmonics. The lowest order term is a radial vector of the form, 



$(x) = 4> a {x)T a 



(r)—T a . 



(4.2) 



2 Each of these regions is bounded by horizons. It is standard practice to refer to the first region and 
horizon as 'the Gribov region' and 'the Gribov horizon', a convention we follow in this paper. 
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on which the Faddeev-Popov operator evaluated at one of our spherical configurations 
operator takes the simple form 

(v A $r = - (V + V - V) - + 2^^V 

\ r r z J r r z 

This expression is positive definite provided / > 0. The Grotowski and Schirmer argu- 
ment now proceeds by evaluating the effect of the Faddeev-Popov operator on all higher 
harmonics in the field, and showing that positive definiteness on radial vectors is sufficient 
to extend the same property to the whole series with only the supplemental constraint 
— 1 < / < 3. Already we see a curiosity - if / can be negative, it no longer seems that ( |4.3| ) 
may be positive definite (see below). This is easily resolved though, as we simply insist 
that < / < 3. With this, positivity of f(r) leads to positivity of T>a which shows that 
the gauge field (and the copy which may be made infinitesimally close to it, as in ( |3.22| )) 
lie within the Gribov horizon. 

4.2 Reinstating the coupling 

We will now add some detail to the above discussion, along with the factors of g. We 
expand a general field 

<I>(x) = ^^ m (r)^ m (0,0)T a , (4.4) 
jlm 

where Y? lm (8, <j)) are an orthonormal basis for su(2) valued-functions on the 2-sphere carry- 
ing orbital angular momentum (j) and total angular momentum (/). The third component 
of the total angular momentum is m. The sum is understood to be constrained by \m\ < I, 
j > and \j — l\ <1 for j > 0, I = 1 for j = 0. The first term in the series has j = m = 0, 
I = 1 and is the radial field (|4.2|) . The Faddeev-Popov operator evaluated at one of our 
spherical configurations and acting on radial vectors gives 

/ 2 2 \x a f — 1 x a 

{V A <$>) a = - Ug 10 + -<f/ 01Q - -3 0oio J — + 2g— —0OW — 

-jd r r d r + 2g 3 — 

The matrix element of the Faddeev-Popov operator between such fields is given by 



— ( - ^d r r 2 d r + 2g- — I <f> 010 . 



J 



d 3 x <t> a {V A (i>) a = 47rjdr r 2 <f> 010 ( - ^d r r 2 d r + 2g f ~ ±± 9 ^ Q1O . (4.6) 
o 

The sign of this expression determines the positivity of the Faddeev-Popov operator. The 
radial Laplacian is a positive definite operator (for suitable boundary conditions), so the 
sign depends crucially on the /-dependent multiplicative term. Since we may in principle 
always find a field for which this matrix element is dominated by contributions from the 
multiplicative term, if / — 1 + g~ l becomes negative the matrix element also becomes 
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negative. For example, take g = 1 and a(r) again as in ( |3,18| ). The corresponding / is 
given in Figure |3] and becomes negative. Now choose </>nio( r ) = ol(t). The field <I>(x) = 
(poio(r)x a /r has finite norm, 



d 3 x (j) a (x)(j) a (x) = 47r / drr 2 a(r) 



47T 

9V3' 



The matrix element of the Faddeev-Popov operator acting on this (j) is 

00 

'2/(r) 1 



(®\V A \$>) = 4vr / drr z a(r) 



^d r r 2 d r )a(r) 



= 4, /d. r» „ (r) (M^M _ /WW)) 


~ -0.85. 

and a plot of the integrand is given in Figure fjj. Thus the Faddeev-Popov operator is not 
positive definite. 

We have now seen that the positivity 
of the lowest term in the harmonic expan- 
sion depends on the sign of the combina- 
tion 



f(r) :=g[f( r )-l+l 
9 



(4.7) 



We now turn to the higher harmonics. 
Writing A? (x) in terms of /, the Faddeev- 
Popov operator becomes 

V A = -V 2 - [A,di-], 

where the gauge field A is defined by 




Figure 7: The density of the matrix element 



AUx) :-- 



Rr) ~ 1 x b 



£iab— = gAi(x). 
r 



All the factors of the coupling have been tidied away into /. The remainder of the Grotowski 
and Schirmer argument now follows with f(r) in [10| replaced by /. So, we divide the matrix 
element ( $ \T>a \ 3>) into four pieces, I, II, III and IV. The first contains the j = term 
which is given above, the remainder correspond respectively to j > 1 with I = j, I = j — 1, 
I = j + 1. The second term is 



m, l=j ' 



1 



II = 4tt / drr2 ^«m(r) ( - ^>, r 2 <>r + 



1(1 + 1) + f(r) - 1 



4>llm{r). 



(4.8) 
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In order for II to be positive we require f(r) > — 1. However we have already seen that 
/ > is required for positivity in the radial part. The third term is 

oo _ 

III = 4vr [ drr " <W r ) ( ~ ^d r r 2 d r + l l±lzIM \ ^( r ) , (4.9) 

m,l=j-l{ ^ T ' ' 

and since j > 2 (for j = 1 the last term vanishes) we require / < 3. The final term is 

00 

IV = 4n Y, j drr " *iim(r)(- ^d r r 2 d r + m ~ 1 + M U fa ( r ) , (4.10) 

m,l=j+l{ ^ r ' ' 

which is non-negative as I > 2 here. Therefore, for < / < 3 each term in the Faddeev- 
Popov matrix element is positive, giving us a positive definite operator. This is the exten- 
sion of the Grotowski and Schirmer argument to arbitrary g. 

4.3 Outside the horizon 

As we have demonstrated, if / becomes negative then we may find fields (we gave a radial 
field as an example) such that the Faddeev-Popov operator fails to be positive definite. 
In this section we will show that all / must somewhere become negative (for any g), and 
therefore the corresponding gauge fields must lie outside the Gribov horizon. 
Let f(r) and a(r) be given by 

= rV(r) + W(r) + 1 _ 
sin(2ga(r)) 
- gr 2 a"(r) + 2rga' (r) 
^ K r > = S m(2ga(r)) ( 4 - n ) 



ga(r) 



titt + Ar as r — > 0, 
mil + /i/r 2 as r —> 00. 



We will now prove that all such f(r) must take negative values. To begin, if ga(r) interpo- 
lates between different multiples of tt at and 00 then at some point it must cross an odd 
multiple of ir/2. We have already seen that in order for f(r) to be well behaved at such a 
point, call it c, then we must have 

7T 

ga(r) = (2n + 1)~ +g-y(r), 
c 3 



7(7") = const, yr ^ j + • • • near r = c. 

At such a point, /(c) is given by 

/(c) = lim r2 gV / ( r ) + 2r 9l'( r ) _ lim r 2 gj\r) + 2rgj(r) 
r^c sin((2n + l)ir + 2g r y(r)) r^c s'm(2g^(r)) 

= -1. 

Thus f(r) must take negative values in a region around r = c. Taking, for example, 
g = c = 1, the resulting /(r) as a(r) passes through ir/2 at r = 1 is given in Figure || 
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We may now restrict to a(r) which 
take the same value at r = and r = do. 
Expand ga(r) = mix + g"f(r) where 7(7*) 
vanishes at r = 0, r = 00. The corre- 
sponding gauge configurations are 



r 2 gj"(r) + 2rgj'(r) 
sin(207(r)) 



(4.12) 



which is independent of m. All such gauge 
transformations lead to the same config- 
urations and copies. If these ga(r) cross 
(multiple times) the value tt/2 then they 
must behave as above and become nega- 
tive. Therefore we need consider only the 
case where a{r) obeys 



a(0) = 0, a(oo) = 0, 



. 9 
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Figure 8: Example of the behaviour of f(r) at 
points (here r = 1) where ga becomes a multiple of 
tt/2 



\ga(r)\ < tt/2. 



(4.13) 



An example of this kind is given by our staple, ( 3.18| ). We know that for such a(r) when 
r is small a{r) ~ ar with the constant a free. When a > there must exist a point r* 
such that a(r*) is a positive maximum; this follows from the behaviour at large r where a 
behaves like b/r 2 and demanding that a be continuous. At the turning point 



/>*) 



r* 2 ga"(r*) 
sm(2ga(r*)) 



(4.14) 



The second derivative is negative as this is a maximum. Since a(r*) > and by assumption 
< 2ga < 7T this implies sin(2ga) > 0. By continuity f(r) is then negative at and in a 
region around r = r* . 

If at small r, a{r) ~ ar with a < then there must exist an r* such that a(r*) is a 
negative minimum. In this case the numerator of ( |4.14| ) is positive but the denominator is 
negative and f(r*) is again negative. 

This completes the proof that the f(r) we construct always take (both positive and) 
negative values. One may imagine that there are other ways to take the limit ( |3.14| ), for 
example by adding higher powers of r or a(r) which vanish in the limit. Such terms can not 
affect the leading order behaviour of a(r) as r — > and 00, however, as these are necessary 
for the limit to exist and f(r) to be well defined. As it is only this leading order behaviour 
which we used to prove that the f(r) become negative, such corrections do not affect our 
argument. 

In summary, we have seen that no spherical copies constructed in this way lie within 
the Gribov horizon, contrary to the statements made by Grotowski and Schirmer, as the 
Faddeev-Popov operator at such configurations fails to be positive definite when acting on 
radial vectors. 
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4.4 Monopole configurations 



As mentioned earlier, it is not physically required that our configurations have finite norm, 
only finite energy. Without the restriction of finite norm we may generate copies of con- 
figurations which behave asymptotically as monopole solutions (we are here thinking of 
embedding our configurations in a Georgi-Glashow type model). We begin with the spher- 
ically symmetric configuration 



AUx) 



1 x° 
gr r 



(4.15) 



which we have written in terms of f(r). We will ask that a(r) and f(r) behave as before 
for small r, but that f(r) — > as r — ► oo so that we asymptotically approach the usual 
monopole profile. This change in boundary condition requires a change in the asymptotic 
behaviour of a(r). It is not difficult to show that if ga{r) = c/r+ga(r), where a(r) behaves 
as 1/r 2 and higher terms as before, that the asymptotic behaviour of f(r) becomes 



r 2 ga" + 2rga' 
sin(2ga) 

sm(2ga) \r 



(4.16) 



which vanishes as r — * oo as required. 

Note that such configurations, like those 
found above, have finite energy. An ex- 
ample is to take 



a(r) 



1 + r 2 



(4.17) 



which at g = 1 has energy \\Fa\\ — 10.6. 
The corresponding gauge configuration is 
shown in Figure ^. All our earlier argu- 
ments extend immediately to this class of 
configurations; the copies may be made 
infinitesimally close, f(r) must take neg- 
ative values and so such gauge fields lie 
outside the Gribov horizon. 




Figure 9: The gauge field f(r) constructed from 
( 4.17 ), with g = 1, which behaves asymptotically 



as a monopole and has Gribov copies. 



5. The loss of colour 

We now return to the construction of colour charges started in Section 2. Recall that 
the ability to associate colour with the charge ( p.9| ) relied on it being gauge invariant. 
This followed from the transformation properties of the dressing ( p. 10 ) which in turn was 
derived from the uniqueness of the Coulomb gauge fixing. The fact that we have copies 
will invalidate this argument. In this section we want to make this breakdown of colour 
explicit. 
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A question not clarified in Section 2 is: what colour does the charge (|2.9| ) have? At 
first sight this looks quite complicated as ty(x) is a composite operator constructed out 
of both the matter and gauge fields which can both contribute to colour. However, the 
overall colour is quite simple to determine in practice. Before we do this, we should say 
that in the bulk of this paper we have described A(x) in terms of a classical field, rather 
than a quantum operator. We are thinking of such objects, and the matter fields, as 
the eigenvalues in the Schrodinger representation. Wave functionals £[^,,-0] = {Ai,ip\ C) 
represent quantum states | () on which Ai(x) and ifj(x) are diagonalised. We satisfy the 
canonical commutation relations by representing the conjugate momenta Ef and ^' by 
functional differentiation. The dressed operator h~ l ijj is therefore also diagonalised. In 
this way we may perform the manipulations below, understanding them to represent the 
action of operators on wave functionals. 

Neglecting copies for the moment, gauge invariance of *$>(x) implies that its colour is 
constant on any gauge orbit. Making this more explicit, we write 

^[A,ip](x) = hr x [A](x)^{x) , (5.1) 

then gauge invariance means 

■$[A U ,ip u ](x) =h- 1 {A u }(x)4> u (x) =h- l [A](x)UU- 1 i){x) = ^[A,tp)(x). (5.2) 

On any orbit there exists a configuration where the Coulomb gauge condition holds and 
hence the static dressing becomes the identity. At this point the colour charge is just 
the matter field and we can easily arrange that it is, say, a localised blue charge there. 
Anywhere else on the orbit the dressing will not be trivial but the overall colour will still 
be blue. In this sense, the dressing is transparent to the matter's colour. As we make a 
gauge transformation away from the Coulomb gauge configuration the matter starts to lose 
its colour but the dressing compensates so that the overall colour is preserved. That is, as 
long as the dressing is sensitive to the applied gauge transformation. 

The existence of copies, though, means that we have gauge transformations for which 
the dressing is insensitive. In particular, the dressing evaluated on the configuration ( p.l| ) 
is unity as the gauge field is already in Coulomb gauge, 

V[A,il)](x) = h~ l [A](x)i){x) =tp(x). (5.3) 



Performing the gauge transformation ( |3.2[ ) parameterised by a(r) obeying ( 3.17| ), the gauge 



field is transformed to its copy ( p.3[ ) , which is also in Coulomb gauge and so the dressing re- 
mains unity. This does not compensate for the gauge rotation of the matter field generated 
by (U|), so that 

x a 

^\A U ,ip u ](x) =U~ 1 ip(x) = cos (ga(r))ip(x) -i sin (ga(r)) — a a ip(x). (5.4) 

r 

It is clear that the commutator of this transformed ^ with the colour charge operator will 
not be the same as that of the original ^, which has therefore lost its colour. 

The copies may be made infinitesimally close together, yet this breakdown of colour 
is non-perturbative, as demonstrated by the coupling dependence of our explicit examples. 
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However, it should be possible to develop a perturbation expansion around such a solution 
and see the breakdown of colour in perturbation theory. As a starting point for this and 
other investigations there are other classes of configurations in the literature where non- 
perturbative effects becomes apparent upon the inclusion of the coupling. We illustrate 
this in Appendix C with Henyey's axially symmetric solutions. 



6. Conclusions 

There are two main conclusions from this paper. The first is that explicit examples of 
Gribov copies in Coulomb gauge are readily produced. These copies arise not only for 
regular gauge transformations that are smoothly connected to either the identity mapping 
or the centre of the gauge group, but also for regular transformations with mixed asymptotic 
conditions which we saw generated open boundaries. Let us re-emphasise that this is a 
wide class of copies - for the case of small gauge transformations any suitably bounded 
function a(r) with the small and large r behaviour 

a(r) ~ r as r — > and a(r) ~ — as r — > oo 



gives a finite norm, finite energy configuration and copy as in ( |3.lD , (|3.13[) . Al- 

though the copies may be infinitesimally close to each other the configurations are non- 
perturbative. This point was only apparent after inclusion of the coupling, which we feel 
it is important to maintain in future studies. 

The other major conclusion from this analysis is the impact these copies have on the 
construction of colour charges and confinement. Although one may define colour pertur- 
batively, such a definition breaks down non-perturbatively. We constructed a matter field 
on which a physical colour was well defined provided that the gauge fixing condition, used 
in the dressing, picked out a unique representative on each gauge orbit. It is here that 
the copies enter, the existence of which contradicts this uniqueness and gives the colour a 
gauge dependence. This shows that coloured charges cannot correspond to physical asymp- 
totic sates and are therefore confined to colourless bound states. This has disentangled the 
fundamental reason for confinement from the detailed dynamics which dictate the scale of 
hadronisation. 

An immediate open question is how we may extend our methods to create copies which 
lie inside the Gribov horizon. An obvious choice is to change the form of the gauge fields, 
as in Henyey's class of copies || where spherical symmetry is replaced by axial symmetry. 
Alternatively, we could change the form of the gauge transformation. 

A deeper question to be addressed is that of the physical significance of the Gribov 
horizon. To what extent does it matter whether a given configuration lies inside or outside 
the horizon? There is at present no physical motivation for working inside the horizon - 
such a choice is essentially arbitrary as we could choose to work in any other Gribov region, 
nor does the choice eliminate the problem of copies. We feel that this work and such future 
investigations as outlined above will help towards understanding these important topics. 

Note added in proof: another method to try to avoid gauge fixing is the gauge invariant 
exact renormalisation group, see e.g. [p 
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A. Baker-Campbell-Hausdorff type formula 

In this appendix we collect together some results associated with the Baker-Campbell- 
Hausdorff series used in Section 2. The basic result which is discussed in many quantum 
mechanics texts is that 

e^e-^ = B + [A, B] + I [A, [A, B}} + I [A, [A, [i,B]]]+- 

= e A B, (A.l) 

where we define A n B by ~A°B = B, AB = [A, B] and A n B = ~A{A n - l B). In addition 
we have the result that if 5 is a derivation then 

e A Se -A = ^Z^_ SA (A2) 
A v ' 

It is helpful to recall how this result is derived as we will need to extend it to the case 
where we have an additional derivation acting on the whole of this equation. 
The proof of flAg) rests on the simple identity that 

L e sA Se -sA = _ e sA, §A ) e -sA = _ e sA gA _ (A.3) 

as 

Integrating this we see that 

e A 5e~ A = - f ds e sA (5A)e- sA = - C ds e sA 6A , (A.4) 
Jo Jo 

from which the result follows. 

We now extend this result and introduce another derivation d and calculate d(e A 5e~ A ). 
Using the Leibnitz property of the derivation, we get three terms, two of which involve the 
derivation acting on e A . Using (|A.2j) and (|A.4|) we find that 



l-e A 



d(e A 5e~ A ) = _ 88 A + 



A A 



5 A . (A.5) 



Expanding in powers of A is straightforward and to quadratic accuracy we find that 

d(e A 5e- A ) = -85 A - \[A, 85 A] - \[8A, 5 A] , (A.6) 
which is used in the main text. 
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B. Open example 



In this appendix we construct an 'open' type gauge transformation which obeys U(0) = 1 
and U{x) —* — 1 at spatial infinity. To simplify the construction we will choose g = 1 and 
look for an a Q which is monotonically increasing between and tt and therefore takes the 
problematic value ir/2 once and only once. We demand that a {r) has continuous second 
derivatives in order for f(r) to be continuous, which implies we require three degrees of 
freedom at each patch. We will construct such an a Q by patching continuous functions 
together, our ansatz for which will be: 



a D (r) 



Are 



- br 



< r < -, 
- - 2' 



7T / 1,1 

- + A(r-^r) -<r<4, 



do da dd 
7r+^- + ^ + ^-4<r<oo. 

J'— j^tt 



(B.l) 



We have chosen to patch the functions together at r = 1/2 and r = 4, choosing a Q (l) = ir/2 
- note that although the functions chosen are not monotonically increasing on their own 
they may be made so on their patches for suitable coefficients. Note the choice of the 
second function, determined by ( |3.26 ). 

At r = 1/2 the equations imposing continuity of a (r), a' (r), a' '(r) may be solved to 

give 



96 



b 
A 



82 _ 
17 ~ 

317T 

"96" 



121 + 7^865)e 41 / 17 -^/ 17 
865 ~ 1.363, 
v/865 ~ 0.052. 



5.492, 



2 

17 V 

7T 

" 96 



(B.2) 



The first patching function has a turning point at 1/6 ~ 0.733 which is outside its range 
of definition and is therefore monotonically increasing. The second patching function is 
always increasing since A > 0. Continuity at r = 4 implies 

5067vr 213tt 



d 2 
d 3 

G?4 



.'865 
32 32 

-984vr + 40vrv / 865 ~ 
1600vr - 64vr\/865 ~ 



~ -117.566, 

604.553, 
-886.861. 



(B.3) 



The derivative of the third patching function has no real roots in its range and is mono- 
tonically increasing there. A plot of a (r) is given in Figure [lO| . 

We now turn to the L 2 — norm of our configurations. f(r) is a continuous function 
and is therefore integrable. For large r (i.e. using the third patching function) f(r) may 
be expanded 

f(r) = 1 + 2^-r- 1 + 0(r- 2 ). 

«2 
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The integrand of ||A|| 2 is therefore 
bound by r~ 2 at large r and A has finite 
norm. The copy has norm \\A U \\ given by 
J3.12 ). The behaviour of the copy is not 
severely modified at infinity; cos(2a Q ) — ► 
1 with r~ 4 corrections which can only im- 
prove convergence of the first term. The 
second and third term go to zero like r -4 
and the integral will converge. The copy 
therefore also has finite norm. 

The energy is also finite as we now 
show. The derivative f'(r) is continuous 
except at the patching points where it has 

only finite jumps (in particular it is smooth at r = 1, a Q = it/2) and is therefore integrable. 
With reference to ( |3.10| ) /(r) behaves as follows at the boundaries, 

(/V) - 1) 2 



Figure 10: a (r) which corresponds to the gauge 
transformations approaching —1 at spatial infinity. 



166V + 0(r 3 



(f{r) - If 



l 6 C-2 + 0(r -3 
at 



fit \ ^3 -2 

a 2 



for r small, 
for r large, 

+ 0(r -3 ) for r large, 



and therefore ||-F|| 2 is also finite. The copy has the same energy. 

Although we have not given a closed form expression for a Q {r) we may make this func- 
tion and therefore f(r) arbitrarily smooth at the patching points by requiring continuity 
of the higher derivatives of a a . This is achieved through the inclusion of corrections to the 
patching functions. Requiring a Q to be monotonically increasing is more restrictive, but 
this may be relaxed anywhere except very near a = ir/2 (or multiples thereof) without 
further complication. 



C. Axially symmetric copies 

We review Henyey's derivation || of the equations for the existence of copies of radially 
symmetric gauge fields, here including the coupling. The gauge field and gauge transfor- 
mation are 

A = ia(r,6)(j> a 3 = —2a(r,6)(f> T 3 , U = cos(ga(r))l + i sm(ga(r))k ■ a, (C.l) 
with k = (cos(</>), sin(</>), 0) and a(r, 9) = rb(r) sin(#). The equation to be solved, c.f. (p^q). 



is 



a(r,9) 



1 



b(r)+r 2 sm 2 (e) [b"{r) + ^b'{r) 



(C.2) 



2grsm(6) f3~ l sin (2gr(3b{r) sin(0)) 

Here (3 is the scaling parameter, a — > (3a which causes copies to coalesce. The boundary 
conditions are that b(r) = 6q + t>2 r2 + • • • with 6q 7^ near the origin and b(r) behaves as 
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1/r 3 at infinity. We take Henyey's example, 



b(r) 



K 



3 • 

\2 



(r 2 + l)\ 

In the limit (3 — > the two copies coalesce to the configuration 

15 rsin(6>) 



(C.3) 



a(r,e) 



2g (r 2 + l) 2 ' 



(C.4) 



The non-perturbative nature of this infinitesimal copy is manifest. As proven in 14] this 
configuration lies on the Gribov horizon. 
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